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Using the Boltzmann weights of classical Statistical Mechanics vertex models we define a new 
class of Tensor Product Ansatzs for 2D quantum lattice systems, characterized by a strong 
anisotropy, which gives rise to stripe like structures. In the case of the six vertex model we 
compute exactly, in the thermodynamic limit, the norm of the ansatz and other observables. 
Employing this ansatz we study the phase diagram of a Hamiltonian given by the sum of XXZ 
Hamiltonians along the legs coupled by an Ising term. Finally, we suggest a connection between 
the six and eight- vertex Anisotropic Tensor Product Ansatzs, and their associated Hamiltonians, 
with the smectic stripe phases recently discussed in the literature. 
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I. INTRODUCTION 

Low dimensional spin systems constitute one of the 
most active areas in Condensed Matter Physics due to 
the experimental findings and the associated theoretical 
activity. These systems are strongly correlated with very 
rich phase diagrams studied by means of a miscellanea of 
analytical and numerical techniques, among which the 
study of simplified variational ansatzs for the ground 
state (GS) and excitations have played a significant role. 

In ID there is a plethora of variational ansatzs: the 
AKLT states [Q, the finitely correlated ansatzs p), the 
Matrix Product Ansatzs (MPA) [g, [| M, the Recurrent 
Variational Ansatz (RVA) ||, etc. All these ansatzs have 
a common structure for the GS wave function which is 
given by the sum, over some auxiliary variables, of prod- 
ucts of amplitudes that also depend on the spin variables 
at the sites. The basic quantity here is the "matrix prod- 
uct amplitude" -A a)l g[mj] where mi is the spin at the i th 
site and a and (3 are auxiliary variables, which can be 
associated to the links meeting at the site. For a spin 
chain with N sites and periodic boundary conditions the 
corresponding state can be written as M 



I^mpa) = y^Tr(^4[mi] • • .A[m N ])\mi) . . 



\m N ) 



(1) 



where the trace is over the auxiliary variables a. Some 
MPA states, such as the AKLT ones, are exact ground 
states of a Hamiltonian, which is given by the sum of 
projectors between nearest neighbours sites Q. In other 
cases the MPA states are used as variational ansatzs for 
a given Hamiltonian, with the MPA amplitudes A a> ^[m] 
playing the role of variational parameters. Within the 
latter category fall the DMRG states (for a review 



on the DMRG see g, |[), which are in fact MPA states 
with open boundary conditions and position dependent 
amplitudes ( i.e. inhomogenous MPA's) j|, [k| O. In 
the DMRG the auxiliary variables label the states kept 
in the blocks. 

The MPA states can be generalized in a natural way 
to 2D systems, replacing the matrix amplitudes A aj( 3[m] 



by "tensor product" amplitudes A ai , 



ml, where z 



is given by the coordination number of the lattice, i.e. 
z = 3 for an hexagonal lattice, z — 4 for a square lattice 
and so on so forth @ fl| |l|, [l§ |l|, [H], [If |l|, |o). 
By analogy with Statistical Mechanics (SM) these states 
can be called Tensor Product Vertex Ansatzs (TPVA) 
because the auxiliary variables on are associated to the 
links of the lattice, while the amplitudes are associated to 
the vertices Another class of 2D ansatzs is formed 
by the Tensor Product Face Ansatzs (TPFA), where the 
amplitudes are associated to the faces of a square lattice 
as in the Face or Interaction Around a Face (IRF) models 
in Statistical Mechanics |l5|, p2] |. The 2D generalizations 
of the AKLT states for spins 3/2, 2 and higher belong 
to the TPVA class. The recipy to construct TPA's is as 
in equation (^), where the contraction of the auxiliary 
variables follows the pattern of the underlying vertex or 
face models. 

Most of the TPA's studied in the literature are 
isotropic, meaning that their properties are largely inde- 
pendent of the spatial direction. However in 2D and 3D 
there are physical systems, like some high-temperature 
superconductors , quantum Hall systems , or man- 
ganites [ p5[ , which exhibit strongly anisotropic properties 
due to the existence of stripes. These objects are static or 
dynamic charge inhomogeneities, which are linear in 2D 
or planar in 3D. One may wonder whether these systems 
can be modelled with simple TPA's, just as Haldane spin 
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chains can be easily described as valence-bond states. In 
this work we shall not address directly this question, but 
the results we have obtained suggest the possibility of a 
simple description of stripes in terms of TPA's. More pre- 
cisely, in this paper we shall investigate a class of TPA's 
based on classical exactly solvable 2D vertex models with 
strong anisotropic properties reminiscent to the stripe 
systems investigated in reference p7fl . Any classi- 
cal SM 2D vertex model, not necessarily integrable, de- 
fined by its Boltzmann weights, give rise to an Anisotopic 
Tensor Product Ansatz ( ATPA) . If, in addition, the SM 
model is exactly solvable, then the corresponding ATPA 
becomes quasi exactly solvable. The latter term is bor- 
rowed from the theory of spectral problems associated to 
the Schrodinger equation psj meaning, in our context, 
that some quantities, as the norm of the ATPA's and 
some expectation values, can be computed exactly in the 
thermodynamic limit. 

To illustrate our proposal we have choosen the well 
known 6 vertex model, whose ID quantum mechanical 
counterpart is the XXZ model or the ID spinless fermion 
plj . We shall show that the corresponding ATPA has 
some similarities with the striped states of 2D spinless 



fermions studied in the literature 1 26 , [27] . 

The organization of this paper is as follows. First of 
all we review briefly the basic ingredients of vertex mod- 
els in Statistical Mechanics ( section II) and the tensor 
product vertex ansatzs ( section III). In section IV we 
introduce the ATPA's based on SM vertex models and 
study their general properties. The ATPA associated to 
the six- vertex model is used in section V as a trial ground 
state for an anisotropic Hamiltonian closely related to 
the XXZ spin chain Hamiltonian, and derive the phase 
diagram. In section VI we briefly comment on the eigh- 
vertex ATPA model. The possible connections between 
the six and eight- vertex ATPA is explored in section VII 
and finally in section VIII we state our conclusions. In 
Appendices A and B we collect some technical results. 



II. VERTEX MODELS IN STATISTICAL 
MECHANICS 







(2) 



The statistical weight of a global configuration is given 
by the product of the Boltzmann weights of all the ver- 
tices. The partition function Z is the sum of these 
weights over all the link configurations, which can also 
be expressed using transfer matrices. The row-to-row 
and column-to-column transfer matrices are defined as, 



a' s i—1 
N 

f's i=l 



(3) 



where £ = . . . , £i), a = (a±, . . . , a^r), etc and the 
periodic boundary conditions are assumed along both di- 
rections. Using (0) the partition function Z reads, 



Z = Tr T, 



N 



(4) 



As an example we display in table 1 the Boltzmann 
weights for the allowed vertex configurations of the six- 
vertex model. The link variables take on two values, say 
and 1, which in the standard notation correspond to 
the right and up pointing ( for a — 0), and left and 
down pointing (for a = 1). The allowed configurations 
satisfy the ice rule a + £ = + r], and the Boltzmann 
weights are invariant under the reversal of all arrows, 
which leaves three independent ones, called a, b and c. 
The 6- vertex model is integrable: there is a uniparamet- 
ric family of transfer matrices T(u) commuting among 
themselves. This is guaranteed by the Yang-Baxter equa- 
tion satisfied by the Boltzmann weights. 



Throughout this paper we shall follow closely Baxter's 
book | pl| . Let us consider a rectangular lattice with N 
rows and L columns. Throughout these paper we shall 
also use the term "legs" for the rows and "rungs" for the 
columns. In a vertex model there is a local state vari- 
able a associated to every link and a Boltzmann weight 
associated to every vertex •, which depends on the four 
link variables meeting at it. We shall represent the Boltz- 
mann weight as 



weight 


a 


b 


c 


V 

a (3 

t 










1 1 






1 
1 


1 

1 1 
1 


1 



1 


1 

1 




Table 1. Boltzmann weights of the six- vertex model 



III. TENSOR PRODUCT VERTEX ANSATZS 



As in the previous section we shall consider a lattice 
with TV legs of length L. In the quantum spin model 
there is a spin degree of freedom m at each vertex of 
the lattice. To construct a TPA we shall associate an 
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auxiliary variable a to each link, as in the SM models. 
The TPA amplitudes will be denoted as 

V 

A%[m] = a - to - (3 (5) 
I 

By analogy with SM we shall define the row-to-row 
and column-to-colum transfer matrix amplitudes 

«m=e n<*:« ■'''>*] ( fi ) 

a' s i—1 
N 

< /3 [m]=^n<;!r i K] 

where m = (mi,...,mi) for yl row while m = 
(mi,...,mjv) for ^4 co1 - Using eqs.(||) the TPA can be 
written in two alternative ways, i.e. 

|V)row = E Tr ( AmW M ■ ■ ■ ^ r ° W [m N ]) Imx, . . . , m N ) rov 

III ' 8 

(7) 

|^>coi = E Tr (A col [ mi ] . . . A col [m L }) \m 1: ..., m L ) col 

m's 

These equations are identical to eq. (jl|) , which implies 
that the TPA can be regarded as a MPA through the 
legs or the rungs. The norm of the TPA is given by, 

<V#) = Tr rZ = Tr T c L ol (8) 

where 

t v7m' = E A vl M^fr M ( 9 ) 

m 

T^, sp > = E ^H^VH 

m 

Thus the computation of the norm (||) amounts to that 
of the partition function of a classical SM vertex model 
where the link variables are twice those of the quantum 
mechanical model. 



IV. ANISOTROPIC TENSOR PRODUCT 
ANSATZS 

A Generic Case 

Let us suppose we are given a vertex model with Boltz- 
mann weights Using them, we shall define a ATPA 

model by the equation: 



A^[m] = 5 m , TI W^ (10) 

where the spin variable at each site, i.e. to, is identified 
with the link variable rj. In the case of the six-vertex 
model the shall adopt the convention that to = cor- 
responds to spin 1/2 and to = 1 to spin —1/2. For the 
six-vertex model the corresponding TPA amplitudes are 
given in table 2. 



Amplitude 


a 


6 


c 


V 

a m f3 
t 




+ 





1 + 1 





+ 1 
1 


1 

1 - 1 
1 


1 

0-0 

1 


1 

1 - 




Table 2. Six- vertex TPA amplitudes 



The choice ([H]) is extremely anisotropic since it treats 
on a very different footing the vertical and horizontal 
directions of the lattice. This is the main reason to con- 
sider both the row-to-row and column-to-column transfer 
matrices, which give rice to complementary descriptions 
of the ansatz. In SM models where the leg and rung 
variables run over different sets, one can obtain two in- 
equivalent ATPA's, not related by a 90 degrees rotation. 
In the rest of the paper we shall suppose that the link 
variables are of the same type in both directions. 

Eq. ( [i7j| ) implies a simple relationship between the row- 
to-row TPA amplitude (||) and the row-to-row transfer 
matrix (j^) of the underlying SM model, namely 

A™%[m] = 6 m , v T^ (11) 

which leads to the following row representation of the 
ATPA state (|), 

I, /A \ * T->row rpiow rpxow m \ 

\f/Tow — / J -L mi m2 ^ m2 m3 • ■ ■ J mN , mi |"11> ■ • • 7 "»N/row 

(12) 

where T™™ m2 is the row-to-row transfer matrix (|J) built 
up from the SM Boltzmann weights defining the 

ATPA ( eq.([To|)). The structure of this state is similar 
to the Kramers- Wannier variational state, first proposed 
for the GS of the Ising model |l3|, [2{|, where the ana- 
logue of T row is played by 2 x 2 matrices. The ATPA 
built from the choice ([Toj ) can be seen as a superposition 
of leg states connected through the row-to-row transfer 
matrix of the SM model. For example, in the antifer- 
roelectric phase of the six-vertex model, the state along 
the legs will be mostly of Neel type and correlated anti- 
ferromagnetically with their nearest neighbours legs. In 
the spinless fermion picture the latter state is a Wigner 
crystal with CDW order. 
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The norm of (|12|) is given simply by 
<V#) row = £(T— m J 2 (T— m3 ) 2 . . . (T— mi ) 2 (13) 

m' s 

It is important to notice that (|l3| ) is not the parti- 
tion function of the SM model defined with Boltzmann 
weights W^'jt or their square. The reason being that in 
general, 

(14) 

where the LHS of this equation is the square of the el- 
ement TI? W ' of the row-to-row transfer matrix, while 

III ] . Ill 2 

the RHS is the entry (mi, m 2 ) of the square of the row- 
to-row transfer matrix. In any case, the computation of 
( |l3| ) requires much less effort than eqs.(^) because the 
matrices involved contain half of the indices of those of 
the general case. In other words, the ATPA does not lead 
to a doubling of indices in the row representation. 

The situation improves even further in the column rep- 
resentation. Using eqs.(||) and ( |io| ) we see that the col- 
umn ATPA amplitudes are given by the product of the 
Boltzmann weights on a column, i.e. 

N 

»=i 

where too = m^. Consequently the column-to-column 
ATPA transfer matrix (^]) becomes 

N 

m i—l 

Let us suppose for a moment that we restrict ourselves 
to the "diagonal" sector of T co1 , which is defined by 
the choices a. = a! and (3 = 0. Then (|l^) becomes 
a column-to-column transfer matrix (|J) with Boltmann 
weights being the square of W^, namely 

N 2 

r^tpp = E II (<f;£-x) = T ^(w 2 ) (i7) 

m i=l 

For a generic TPA this diagonal truncation may be a 
good approximation in certain regions of the parameter 
space, as has been shown by Niggemann et al. in a TPA 
for a spin 3/2 system on a hexagonal latice [O]. 

We shall show below that for a subclass of ATPA's, this 
diagonal truncation is in fact exact, which has important 
consequences. 



B Ansatzs with conserved quantum numbers 

Let us assume that the Boltmann weights satisfy 
a conservation law of the type, 

W a,(! = unlcss a + £ = I 3 + V ( 18 ) 

where the link variables label the basis of an irreducible 
representation (irrep) of a Lie group Q. The six- vertex 
model corresponds to the spin 1/2 irrep of the group 
Q = SU(2), with the convention a — (resp. 1) for 
the s z = 1/2 ( resp. s z = —1/2). For a general Lie 
group the link variables will be given by the weights of 
the corresponding irrep. 

The immediate consequence of ( |l8| ) is that the non 
vanishing terms of ( |l6| ) must satisfy 

a, t + rrii-i = fa + to, ; , i=l,...,N (19) 
+ TOj_i = fa + rrii 

which implies 

a'. - Ul =p'.-fa = Q i (20) 

where Qi = 1,0,-1 for the six- vertex model. In the 
general case Qi, being the difference of two weights of 
irreps, is either zero or a root of the Lie group Q . 
Defining the Boltzmann weights W® as 

{W Q f a } = W% W»S* (21) 

we see from (|l6]) and ( |20| ) that T co1 breaks into block 
transfer matrices J 1 *** labelled by the vector Q = 
(Qi, ■ ■ ■ , Qn), whose entries are given by 

N 

m i— 1 

The case Q = corresponds to the matrix ([l7]), and 
hence the truncation of the model to the "diagonal" sec- 
tor is not an approximation, but an exact result. This 
fact greatly simplifies the computation of the norm of 
the ATPA in the thermodynamic limit L — > oo, which 
is given by A^ ax where A max is the biggest of all largest 
eigenvalues of the matrices T®. In appendix A we 
show that this eigenvalue belongs to the Q = sector 
and thus, 

lim (^} col = A^ iax , A max = A Q=0 (23) 

L — >oo 

This is quite a useful result for it implies that if the 
SM model defined by the Boltzmann weights W° = (W) 2 
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( pl| ) is integrable then A max can be computed exactly, at 
least in the limit N — > oo. In the case of the ATPA based 
on the six-vertex model, the Boltzmann weights W° are 
simply the square of the original ones, i.e. 



W°(a,b,c) = W{a 2 ,b 2 ,c 2 ) 



(24) 



and hence the norm of the ATPA can be computed ex- 
actly in the thermodynamic limit. 

There is yet another important consequence of the con- 
servation law (|l8|). As it is well known in the theory 
of transfer matrices, eq.(|l8|) implies that the row-to-row 
transfer matrix preserves the sum of all quantum num- 
bers of every row, i.e. 



I^M+ll^row = E( T mTm 2 ) 2 (^ 2 W m a) 2 



(27) 



//-prow \ 2 /Ad 

■ ' III N HI 1 / 



It can be shown that the square of the row-to-row 
transfer matrix can be written as ( recall eqs.(^p|) and 

©) 



row \2 
mi ,rri2 / 



Q 



(28) 



and hence the sum (|27| ) becomes 



L L 

T m,m' = UllleSS m i = m 'i ^ 

1=1 i=l 

Hence all the terms appearing in the sum Jl^), giving 
|^)row> must have the same value of "angular momenta" 
per leg. In the six-vertex model this implies the vanish- 
ing of all correlators between raising and lowering spin 
operators among different rows/legs, i.e. 

(St,a S7, b ) =0, if 3 (26) 

where S i is the raising (lowering) spin operator on the 
a th site of the i th leg. In other words, the quantum fluc- 
tuations across the legs of the ATPA are strictly forbid- 
den. In the six vertex model the spins may only fluctuate 
along the legs. Using the spinless fermion terminology, 
the only allowed charge fluctuations occur inside the legs. 
As mentioned in the introduction, this lack of quantum 
fluctuations across the legs is reminiscent to that occur- 
ring in some models of high-T c superconductors ( see sec- 
tion VII). 

The previous considerations give us a hint on what 
sort of Hamiltonians the ATPA's may become approxi- 
mate ground states. After all, we want to use the ATPA 
as variational ansatzs for physically interesting systems. 
We postpone this question until next section after a dis- 
cussion on correlators and density matrices for ATPA's. 

C Correlators and density matrices 

Let Of i+l be a diagonal operator acting between the 
legs i and i + 1 that do not change their states and with 
matrix element C^. m , +1 ■ A typical example in the six- 
vertex model is provided by the Ising term cr* a af +1 a . 
The expectation value of Of i+1 in the ATPA is be given 
by, 



<^ +1 |^} ro w = E ]T T(W^) mi , m2 

m'sQi,...,Q N (29) 

In the thermodynamic limit this sum will be dominated 
by the term Q\ = • • • = Qjy = 0, just as in the computa- 
tion of the norm of the state and hence the expectation 
value of Of i+1 reduces to an expectation value in the SM 
model with Boltzmann weights W°. This is a property 
of all diagonal operators which allow their exact evalua- 
tion, provided they are known in the underlying exactly 
solved model. 

In the case of the operator of a of +1 a , its correlator is 
equivalent to the SM expectation value 

P = ( a i,a °f+l,a)ATPA = ( s i,a S i+ i ia ) 6 _ vor tcx (30) 

where s = 1, — 1 is related to the link variable a — 0, 1 by 
the equation s = 1 — 2a. In appendix B we shall compute 
this quantity using the exact solution of the six-vertex 
model. 

Let us next consider an off-diagonal operator 0° d 
acting on the i th leg with matrix elements O 1 ^. , — 
(rrii|0° d |m|). Its expectation value will be given by, 

(Of) = Tr (pi O od ) (31) 

where pi is the density matrix of the i th leg whose entries 
are, 

P mi ,m! = mmT X! ( T m*m 2 ) 2 '-- ( 32 ) 

m's/rrii or m( 

frirOW rpTOW rpTOW rpTOW /rrirOW \ 2 

mi_i ,mi m;_ i ,mj mi ,mi+i mprrij_j_i * ' ' V mjsj ,mi / 

Using again eq.(|28|) in the thermodynamic limit we can 
write pi as 



6 



A 2 

max 



E ■ 

mi_i,m i+ i 



(33) 



ry-irOW rpTOW T 

1 mi,mi+i 1 m!,mj+i ^mi+i 



l/r 

where «m are the left and right eigenvectors with highest 
eigenvalue A max of the transfer matrix T(W°). Eq.(|3^) 
shows that the regions located above or below of a given 
leg behave as if they were in a single coherent state, which 
in some cases can be identified with the ground state 
of the underlying quantum mechanical model. This is 
indeed the case if we assume that T(W°) is a symmetric 
matrix, which is achieved in the six-vertex model if the 
Boltzmann weights a and b are equal. From now on we 
shall assume the latter condition which implies that v l m = 
v r m = v m . 

The computation of (^Tj) is in general quite difficult de- 
pending on the operator in question. An approximation 
can however be made using the following result. If 0° d 
is a positive definite operator then 



{Of) < (v\O oa \v) 



(34) 



The proof of (34) uses the Perron- Frobenius theorem 



and the Schwarz inequality and it is similar to the one 
given in Appendix A to prove eq.(|23|). For diagonal op- 
erators acting on a leg, eq.(M) becomes an equality. 



V. THE SIX- VERTEX ATPA 

The considerations made above suggest that an ATPA 
based on the six- vertex model should be a reasonable ap- 
proximation to the ground state of the following Hamil- 
tonian: 



H — Hleg + -ffrung (35) 
N L 

#lcg = - g EE KXa+l + <a<a+l + A 0< o <a+l) 
i— 1 a— 1 
^ N L 

H lung = — J a i,a a i+l,a 
i—1 a— 1 

which is a combination of the XXZ Hamiltonian along the 
legs and an Ising one along the rungs. The latter choice is 
motivated by the absence of quantum fluctuations across 
the legs. This model has also been studied in reference 
p2| using bosonization techniques. Using the Hellman- 
Feynman theorem and eqs.(|30|) and (|34]), one can find the 
following lower bound of the energy per site of the ATPA 



£(A,A ,J') = £o(A) + (Ao-A) 



dE 
OA 



-J' P(A) 



(36) 



where A is the anisotropy parameter associated to the 
Boltzmann weights W°, i.e. 



A = 



„2 , U2 _ „2 
"0 + °0 c 



2a 6 



(37) 



a = a ,b 



b\c 



P(A) is the expectation value defined in eq.(|30|), and 
Eo(A) is the GS energy per site of the XXZ model with 
anisotropy A. 

The problem is: fixing A and J', find the value of A 
that minimizes the total energy (pq), i.e. 



(38) 



A = A(A ,J') 



It is easy to see that if J' = then A = Ao- 

In the antiferromanetic region (AF) of the XXZ model, 
i.e. A < —1, the parametrization of the Boltzmann 
weights is given by pl| , 



A = -coshA, A > 
ao = p sinn 

bo = p sinh 



2 

\ + v 



(39) 



Co = p sinhA 

where p is an overall factor and v is the spectral pa- 
rameter that is set to zero in order to have a symmetric 
transfer matrix. 

The total energy and its derivative in the AF region 
can be found from the Bethe ansatz solution and they 
read 



1 °° 1 

E a = -coshA - sinhA - 4 sinhA -^-^ (40) 

u 2 ^ e 2mA + 1 

m— 1 

OEo 1,1,4- 1 (41) 



OA 2 tanhA tanhA 



v m e 2mA 



E 



The matrix element P in this region is derived in Ap- 
pendix B and it reads, 



P = -1- 



sinhmA 



sinh 2 (A/2) tanh(A/2) ^ cosh 2 mA 



(42) 



In the critical region (C), i.e. —1 < A < 1, the 
Darametrization of the Boltzmann weights is given by 
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AF / 


b / 


a /y r 


*d 



-6 -5 -4-3-2-10123 

FIG. 1: Phase diagram of the Hamiltonian (|35| ) obtained 
using the six-vertex ATPA. The point a corresponds to 
A = —4.6 and J' = —2.17, b is a generic point on the 
line AF/C above the point a, and finally d corresponds 
to A = 1 and J' = 0. 



-COS/X, < fl < 7T 



u — v 

a = p sin— 
Dq = p sin 



(43) 



Co = p sm/i 

The GS energy per site and its derivative read [Ell 



1 sin/i 
£/ = -cos// 

2 /i 



2 log2 - 2?r , — 

o sinh nx 



dx 



log cosh/xx 



(44) 



<9A 



1 2tt f°° x tanh/ix / 1 
-- A I dx 



A* Jo 



sinh 7rx 



1 



/i tan/i /i z 



(45) 



-21082 + 2^/ bg C0 2 SV " 
/o sinh 7rx 



and P is given by ( see Appendix B), 



P= 1 - 



4~ 



/x tan(/i/2) 7 



<ix 



sinh px 



cosh/ix sinh 7rx 



(46) 



Figure 1 shows the phase diagram of the Hamiltonian 
( |35| ) obtained by minimization of the energy (^). We 
recall that ( p6| ) is a lower bound of the energy of the 
ATPA, and hence does not yield an upper bound of the 
exact GS energy of (|35|). 

The region denoted AF in fig. 1 corresponds to the 
cases where A(Ao, J') lies inside the antiferromagnetic 
regime A < — 1. The region C describes the critical 



regime, i.e. — 1 < A < 1, while the region F denotes the 
cases where A = 1. The phase boundaries between these 
regions have different properties. The AF/C boundary 
line ab corresponds to A = — 1, and hence the transition 
between the AF and C phases seems to be continuous. 
Below the point a the value of A, near the line AF/F but 
on the AF side, is smaller than —1, indicating that the 
AF/F boundary is discontinuous. The C/F line ad is 
also discontinuous, meaning that A jumps across it. Fi- 
nally, there is no discontinuities across the C/F boundary 
above the point d. 

In reference |32| the Hamiltonian (^) was studied us- 
ing bosonization, mean field and renormalization group 
(RG) methods. Disregarding the inter- leg forward scat- 
tering and umklapp terms, that arise upon bosonization, 
the main conclusion of reference J3^| is the existence of 
an AF region whenever | J' | > 2 Ao (mean field result) or 
J' | > 4Ao ( RG result). Furthermore, the inter-chain 
forward scattering terms can be taken into account |33]| 
usingthc sliding Luttingcr liquid approach of references 
p6| , \F\ . In |33| it was shown that the effect of the inter- 
chain forward scattering is to modify the phase bound- 
aries separating the AF and C regions in an asymmetric 
way. Indeed the AF region appears when J' > C+Ao ( 
if J' > 0) and - J' > C_Ao ( if J' < 0), with C+ ^ C-. 



Hence the results of references ]32j, 33 , which should be 
valid in the weak coupling regime A |, | J'\ « 1, suggest 
that the system should be in an AF phase whenever the 
legs are antiferromagnetic, i.e. Ao < 0. This is in contra- 
diction with the ATPA result where there exist critical 
region with Ao < 0. On the other hand the ATPA agrees 
with the aforementioned works on the existence of large 
regions in the phase diagram where the system is critical, 
which we identify with the sliding or smectic Luttinger 
liquid fixed points of p6], pTj] . 



VI. THE EIGHT- VERTEX ATPA 

An ATPA closely related to the six-vertex one can be 
built from the Baxter's eight-vertex model, whose Boltz- 
mann weights are those of the six-vertex model plus two 
new weights Wqq = VF™ = d [|l] . The conservation law 
(plf) now becomes, 



W, 



0,v 



unless a + £ = f3 + r/, (mod 2) (47) 



The transfer matrix T c ° l also breaks into block matri- 
ces with the difference that Qi only takes two values 
and 1, since Qi = — 1 = l(mod 2). 

The eight-vertex ATPA can be taken as the an ansazt 
for the GS of the following Hamiltonian: 
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H — Hi eg + H IUDg (48) 

H\cg = ~2 XI { J x°i,a°i,a+1 + Jy^ta^ta+l + J *°f,a°ta+l) 
i,a 

-ffrung = - J a i,a°i+l,a 



The phase diagram of this model can be worked out us- 
ing the Baxter's exact solution of the eight- vertex model, 
as we did for the six-vertex one in the previous section. 
The results will be presented elsewhere. 



VII. THE ATPA AND STRIPES 

An interesting feature of the six and the eight-vertex 
ATPA's is their possible connection with the stripes in 
high T c superconductors, specially when regarded as elec- 
tronic liquid crystals |^(| g7], [SO) . In this section we shall 
briefly explore this issue which deserves a more detailed 
study in the future. 

The first observation is that the eight-vertex Hamilto- 
nian ([l8]) ( and similarly the six-vertex one ( |35| ) ) can be 
Jordan- Wigner transformed onto the following spinless 
fermion Hamiltonian, 



h = — Yl [ + J »)(i^i.» + h - c -) 

i,a 

+ (Jy - Jx){^i,a^i+l,a + h.C.) 

1 1 (49) 

+ 2J z (n ha - -)(n l+1:a - -) 

-2J'(n 4 , a - -)(n i;a+1 - -) 

which describes the motion of holons along the legs of 
a 2D lattice ( term (J x + J y )), which are coupled by 
density-density interactions ( term J'), together with pair 
tunneling between the legs and the environment ( term 
(J x — J y )) Upon bosonization eq.(49) has a struc- 

ture similar, but not identical, to the "smectic" Hamil- 
tonian in the spin gap case considered in reference |2q| 
and the spinless sliding Luttinger model of reference [B7) . 
Indeed, the smectic symmetry <p a — > <\> a + a a p6[ , where 
4> a is the boson field of the a th leg, is the dual version 
of the standard U(\) symmetry of the six-vertex model, 
which corresponds to 6 a — » 6 a + a a , where 9 a is the dual 
boson fS. The CDW coupling among the stripes in 
pq , pTf corresponds to the term J', while the Josephson 
tunneling is somehow reflected by the pair creation and 
annihilation terms. Assuming these correspondences it 
is quite natural to conjecture a relationship between the 
smectic phases of references |2(| |2?| and the correspond- 
ing phases of the eight-vertex model. The stripe crystal 
phase should correspond to the antiferromagnetic phase, 



the smectic superconducting phase should correspond to 
the disordered phase and finally, the smectic metal should 
be associated to the critical phase, which is the one of the 
six- vertex model when — 1 < A < 1. 



VIII. CONCLUSIONS 

In this paper we have proposed a new class of 
Anisotropic Tensor Product Ansatzs (ATPA) using the 
Boltzmann weights of classical Statistical Mechanics ver- 
tex models. 

We have shown that the computation of the norm and 
some observables simplifies enormously, becoming exact 
whenever the underlying SM model is exactly solvable. 

The strong anisotropy of the ATPA's is reflected in the 
absence of quantum fluctuations across the legs of the 2D 
lattice, a property which suggests a possible connection 
with some current models of stripes. 

We have studied the ATPA based on the six-vertex 
model, as a trial state for the ground state of a Hamil- 
tonian given by the sum of XXZ Hamiltonians along the 
legs of a 2D lattice, which are coupled by an Ising term. 
Using the exact solution of the six- vertex model we have 
proposed the phase diagram of this model and compared 
it with the one obtain with other methods [[32], |33| . 

We have suggested a connection between the six- vertex 
and eight-vertex ATPA's, and their associated 2D Hamil- 
tonians, with the smectic stripe phases considered in ref- 
erences |26|, |27j . 

Let us finally comment on the relation between the 
ATPA and the DMRG. As we explained in section IV, 
the link variables along the legs and the rungs of the 
SM vertex model can be of different type. For example 
we can choose the rung variables £,77 to take only two 
values, say and 1, as in the six vertex model, while 
the legs variables a, (3 can take a large number of val- 
ues, say 1,2,..., 771, as in the DMRG. The ATPA so con- 
structed would have a spin 1/2 at each site with strong 
correlations along the legs. This state would be a sort of 
anisotropic DMRG state with a stripe like structure built 
in. The problem is to device an algorithm to update the 
local weights. 
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APPENDIX A: HIGHEST EIGENVALUE OF T Q 



In this appendix we shall give a proof of eq. (Eq) under 
the condition that all the TPA amplitudes are non nega- 
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tive. Let us call the largest eigenvalue of the transfer 
matrix defined in eq.(22). The statement is that 



APPENDIX B: THE TWO POINT CORRELATOR 

P(A) 



A?<A^ =0 , VQ (50) 
with positive entries and 



Choosing two vectors x an d 
scalar product equal to 1, 



(x\4>) = i 

Xa=U 2 a , (f> a =V 2 a , U a ,V a >0 

one has by the definition of Aq , 



(51) 



(x|r Q |^> <a 



Q 



(52) 



where the equality holds whenever x an d 4> are the left 
and right eigenvectors of respectively (recall that, 
by the Perron- Frobenius theorem, the eigenvector of 
with highest eigenvalue, has all its entries positive). 
Using (B) and (02) we can write the LHS of (j52[) as 



In this appendix we indicate how to compute the ex- 
pectation value 



P (Si,a ^2+1, a)6— vertex 



(58) 



in the six- vertex model with Boltzmann weights a, b and c 
when a = b, which is the case under study. This quantity 
is similar, but not identical, to the polarizability Pq = 
(ai) defined by Baxter 

The partition function of the six-vertex model can be 
expanded as 



Z 



'y ' a "l+»2 foria+rn c n 5 +n B 



(59) 



where n\ and rii are the number of vertices with Boltz- 
mann weight a, etc. The weights a and b contribute to 
P with +1 while c does it with —1, hence P is given by 
the formula, 



{ X \T^)=J2^>? (53) 
= uivlA^ p [rn]A^ w+Q [rn] 



cn,/3,m 



where / denotes the triple (a,/3,m) and xf stands for 



it turns out that 



On the other hand, the Schwarz inequality 



(54) 



(55) 



lim 

N,L^oc NL Z 



d 



d 



da db 



d_ 
da 
d_ 
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db dc 



Z 



(60) 



where / is the free energy per site in the units ksT ' = 1. 
It is important to realize that the derivatives in ( |60| ) 
are performed keeping the remaining ones unchanged. 
Eq.(|60[) assumes that a,b and c are independent quanti- 
ties, however if a = b the formula for P becomes, 



(61) 



In the AF region the free energy / is given, in the 
parametrization by pl|, 



Qi 



< 



\ 



implies 



E 



(56) 



/ = log a 



X + v 



E 



sinh m(A + v) 
n cosh mX 



(62) 



while in the critical region one has, in the parametrization 



\(x\T Q m<A° Q 



(57) 



Hence, choosing x and <fi the left and right eigenvectors 
of one derives the desired result (I5C 



/ = log a 



dx sinh (A + v)x sinh (7r — fi)x 



sinh ttx cosh/ix 



(63) 



Using eqs.(|T]), ( |62| ) , (|6^), ( ^9|) and (B^), with v set equal 
to zero, one can derive eqs.(|42|) and ([16[) yielding P(A). 

In figure 2 we plot P(A) in the AF and C regions. In 
the AF region one has —1 < P < —1/3, while in the C 
region — 1/3 < P < 1. At the isotropic point A = — 1 
one finds P = —1/3, while in the XY model, i.e. A = 0, 
the result is P = 0. In all the ferromagnetic (F) region, 
i.e. A > 1, one has P = 1. 
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[r 



FIG. 2: Plot of P{A) = (a i>a s i+1)0 ) 6 _ vertex in the AF 
region (eq.(B5)) and critical region (eq.(E6l)). 
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